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Formation of stationary localized states in one-dimensional chain as well as in a Cayley tree due 
to a linear impurity and a nonlinear impurity is studied. Furthermore, a one-dimensional chain 
with linear and nonlinear site energies at the alternate sites is studied and rich phase diagrams of 
SL states are obtained for all systems we considered. The results are compared with those of the 
linear and nonlinear systems. 



I. INTRODUCTION 

The discrete nonlinear Schrodinger 

equation (DNLSE) describes a number of phenomena 
in condensed matter physics, nonlinear optics and other 
fields of physics Jl|-Q. The DNLSE in one dimension is 
written as 



.dC„ 
'~dT 



e n C n + V(C n+1 + C n _i) + Xnfn(\C n \)C n , (1) 



where n = 1, 2, 3, • • • , N. The e„ is the linear site 
energy at the site n, V is the nearest-neighbor hopping 
matrix element, and Xn is the nonlinear parameter as- 
sociated with the nth site of the lattice. Since J2 n l^™| 2 
can be made unity by choosing appropriate initial con- 
ditions, |C„| 2 may be considered to be the probability 
of finding a particle at the nth site in the lattice. Con- 
sequently, C n may be interpreted as the probability am- 
plitude of a particle at the nth site. In the context of the 
condensed matter physics, the nonlinear term f n (\C n \) 
arises from the coupling of vibrations of masses at the 
lattice point to the motion of a quasiparticle in the limit 
of high frequency lattice vibrations. If the lattice oscil- 
lators are purely harmonic in nature and are indepen- 
dent of one another, it follows that / n (|C n |) = |C n | 2 
As an application of DNLSE particularly in con- 
densed matter physics, we cite among others the ex- 
citon propagation in Holstein molecular chain |Q. In 
general, the exciton propagation in quasi one dimension 
having short range electron-phonon interaction can be 
adequately modeled by the DNLSE Other exam- 
ples include nonlinear optical responses in superlattices 
formed by dielectric or magnetic slabs |^] and the mean 
field theory of aperiodic array of twinning planes in the 
high T c superconductors jl0| . 

One important feature of the DNLSE is that it can 
produce stationary localized (SL) states. The SL states 
might play an important role in the nonlinear DNA dy- 
namics |ll| and also in energy localization in nonlin- 
ear lattices jl2| |. It has been shown that the presence 
of impurities can produce SL states in one, two and 
three dimensions jl3| The formation of SL states 
in one-dimensional chain with a single nonlinear impu- 
rity and that with a dimcric nonlinear impurity have 
been studied in detail [|l3|,^0). The same problem has 
also been studied starting from an appropriate Hamil- 
tonian |H]j2l|] . The fixed point of the Hamiltonian which 



generates the appropriate DNLSE can also produce the 
correct equation, giving the formation of SL states. The 
formation of SL states in a perfect nonlinear chain con- 
taining a single nonlinear impurity and also in a per- 
fect nonlinear chain with a dimeric nonlinear impurity 
has been studied 22 1. The formation of SL states in a 
perfectly nonlinear Cayley tree due to various nonlinear 
impurities has also been considered earlier pp| , p3| , p4[ . 

On the other hand, there has been a detailed study 
on the formation of localized states due to linear impu- 
rities in one-dimensional chain as well as in the Cayley 
tree p5| . The number of localized states do not ex- 
ceed the number of impurities in case of linear systems. 
We note that the stationary localized states have ap- 
peared due to the presence of either linear impurities 
or nonlinear impurities. Thus, one may naturally ask 
the question of what happens when both the linear and 
nonlinear impurities are present at the same time in the 
system. In order to investigate the mixed role of linear 
and nonlinear impurities in the formation of SL states, 
we consider a simple one-dimensional chain as well as 
a Cayley tree consisting of a single linear impurity and 
a single nonlinear impurity. Furthermore, we consider 
a one-dimensional chain with linear and nonlinear site 
energies at alternate sites of the lattice. 

This paper is organized as follows. The one- 
dimensional chain with a linear impurity and a non- 
linear impurity is considered in Sec. II. In Sec. Ill, we 
consider the one-dimensional chain with linear and non- 
linear site energies at alternate sites. The Cayley tree 
with a linear impurity and a nonlinear impurity is con- 
sidered in Sec. IV. We finally summarize our observa- 
tions in Sec. V. 



II. ONE-DIMENSIONAL CHAIN WITH LINEAR 
AND NONLINEAR IMPURITIES 

Consider a one-dimensional chain consisting of a lin- 
ear impurity of strength e placed at the zeroth site and 
a nonlinear impurity of strength \ at the first site. The 
relevant Hamiltonian for such a system can be written 

as 
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where the nearest-neighbor hopping matrix element is 
assumed to be unity. As a possible solution for the sta- 
tionary localized states, we consider 



(3) 



where E is the stationary state energy. It is also well 
known that impurity states in one-dimensional system 
are exponentially localized. Therefore, the presence of 
two impurities at the two consecutive sites (site and 
site 1) in the chain, suggests us to consider the following 
forms of <b n : 



K = <j>i [sgn(E)r]} r 
0_i n i = 0o [aga(E)rj\ 



\n\ 



n > 1 
n < 0. 



(4) 



The form of <j> n given above are exact and may be de- 
rived from the Green's function analysis |2(J ■ Here rj lies 
between and 1 and is given as 
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Since we are dealing with the localized states, (f>o and <fii 
can be assumed to be real without forsaking the math- 
ematical rigor. Three possible cases may be considered, 
i.e., 4>i = (f>Q, 4>\ — -~4>o, and 0i ^ tfio, the solutions for 
which are called, respectively, symmetric, antisymmet- 
ric, and asymmetric solutions. To take into account of 
all three possibilities, we introduce a variable (3 = <f>o/<pi. 
Consequently, (3 = 1 corresponds to the symmetric so- 
lutions, f3 = —1 to the antisymmetric solutions, and 
f3 7^ ±1 to the asymmetric solutions. Substituting the 
ansatz of Eq. (4) and the definition of (3 to the normal- 
ization condition ^2 |C«| 2 = lj we obtain 
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Using Eqs. (3), (4) and (6) in Eq. (2), we obtain the 
effective Hamiltonian for the reduced dynamical system: 



(i-v 2 ) 
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We first consider the case of (3 - 
Hamiltonian for this case becomes 



±1. The effective 



H±=sgn(E)r,± 



(l-rf) , e(l-77 2 ) 
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where "+" refers to the symmetric states and "— " cor- 
responds to the anti-symmetric states. For a particular 
system, e and x are constants, while 77 is a variable which 
determines the energy of the localized state. The value 



of 77 should be determined self-consistently and, for this 
reason, we treat 77 as the dynamic variable in the ef- 
fective Hamiltonian of the reduced dynamical system. 
The fixed-point solutions corresponding to the localized 
states can be found, from the variabtional principle, by 
the condition dH^/dr) = with t] e [0, 1] . Since the 
form of the probability profile is rigorous fl25j , we obtain 
the correct values of 77, satisfying the equation 



X 
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(9) 



^From Eq. (9) it is clear that the states will appear 
above the band for > an d below the band for 
X* 1 < 0. [Note that the band edges are at E/V = ±2.] 
We will consider positive values of \^ and, hence, the 
sgn(-E) may be taken to be positive. For symmetric 
states, the value of f+(rj) exhibits the extreme values: 
f+(r]) — > as 77 — > and /+(??) — > 00 as 77 — > l/(e/2 + l). 
Thus, /+(7?) increases as 77 increases from 77 = and 
diverges at 77 = l/(e/2 + 1), and it becomes negative 
for 77 > l/(e/2 + 1). This means that there will al- 
ways be only one solution (symmetric state) for all pos- 
itive values of e and x- F° r antisymmetric states, we 
found that f-{rj) — * as 77 — * and f-{rj) — > 00 as 
77 — * l/(e/2 — 1). The divergence of f-{rf) is, however, 
confined between and 1 if and only if e > 4. There- 
fore, for e > 4, /_ (77) increases as 77 increases from 7/ = 0, 
diverges at 77 = l/(e/2 — 1), and finally becomes nega- 
tive for 77 > l/(e/2 — 1). This implies that there will 
always be one antisymmetric state for e > 4 for any 
value of x > 0. However, the situation is different for 
e < 4; the value of f-(rj) approaches as 77 — * and 
77 — ► 1 and remains positive and finite for < 77 < 1. 
This implies that f-{rj) has at least one maximum in 
the range of < 77 < 1 . We confirm by a graphical anal- 
ysis that there is only one maximum. Therefore, there 
will be a critical value of x separating two antisymmet- 
ric states from the region of no antisymmetric state for 
e < 4. The solid line in Fig. 1 is such a critical line 
and is denoted by xtr- Thus, there are three regions in 
the entire plane of e versus Xi having different number 
of symmetric and antisymmetric SL states. The region 
of e > 4 has one symmetric and one antisymmetric SL 
states, while the region of e < 4 is divided into two re- 
gions by the critical line x%-> below of which has only 
one symmetric state and above of which has one sym- 
metric state and two antisymmetric states. Therefore, 
the maximum number of SL states contributed by the 
symmetric and anti-symmetric solutions is three. 

Next, we consider the asymmetric solutions, i.e., for 
[3 7^ ±1. The effective Hamiltonian in Eq. (6) should be 
considered as a function of two variables, 77 and (3, while 
e and x remain constant for a particular system. The 
fixed-point solutions corresponding to the SL states can 
be obtained from the following two extreme conditions 



dH eS /d V = 
dH eS /d/3 = 0. 



(10) 



After a trite algebra, the fixed-point equations reduce 
to 



2 



(1 + /3 2 )(1 + I3 2 - 2 V I3 - en) - Xf?/3 4 (1 - V 2 ) = 

(1 - /3 4 ) - e/3(l + (?) + x /3 3 (l - v 2 ) = 0. (11) 

These coupled equations can not be solved analytically 
and, accordingly, we analyze them numerically to find 
the number of asymmetric SL states (number of solu- 
tions with r] S [0, 1]). 
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FIG. 1. The phase diagram for SL states in the plane of 
e versus x f° r a system of a one-dimensional chain with a 
linear impurity and a nonlinear impurity. Different regions 
are marked with the number of states and the kind of states; 
the region marked by "ls+2a+2as" means that the region 
has one symmetric, two antisymmetric and two asymmetric 
SL states. 



Fig. 1. The energy of one of the two new states increases 
and that of the rest decreases and approaches to the 
band edge as x increases. Therefore, for x < 10.39, the 
localization length of the state, i.e., the characteristic 
length of the profile, decreases as x increases, whereas, 
for x > 10.39, the localization length of one of the three 
states increases as x increases. Thus, for very large x> 
two of the states are strongly localized and the rest is 
weakly localized. 
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FIG. 2. Variation of energy as a function of x is shown for 
e = 2. 



III. LINEAR AND NONLINEAR SITE 
ENERGIES AT ALTERNATE SITES 



It is found that there are three critical lines separating 
different number of states in the plane of e vs x denoted 
by the dashed curves in Fig. 1. Maximum number of 
asymmetric SL states obtained is two. Considering all 
possible states we find that maximum number of SL 
states is five. The "s" , "a" and "as" stand for symmet- 
ric, antisymmetric and asymmetric states, respectively. 
For example, the region marked as "ls+2a+2as" has 
one symmetric SL state, two antisymmetric SL states, 
and two asymmetric SL states. It should be noted that 
the number of states changes one by one as we go from 
one region to another in the phase space. As an exam- 
ple, as we cross from the region "ls+las" to the region 
"ls+2a+las" , the number of states increases by two; 
however, the number of states on the critical line which 
separates two regions is "ls+la+las", yielding the in- 
crement of the number of states one by one. The maxi- 
mum number of SL states with two nonlinear impurities 
was also found to be five ^3|. We thus conclude that 
one linear impurity and one nonlinear impurity play a 
similar role to the two nonlinear impurities as far as the 
maximum number of SL states are concerned. However, 
the phase diagrams for two cases are different. 

Figure 2 shows the variation of energy of the states 
as a function of x f° r a fixed value of e = 2. As x 
increases, the energy of the state increases (goes away 
from the band edge) and, beyond x ~ 10.39, two more 
states appear as expected from the phase diagram in 



In this section, we consider a one-dimensional chain 
with alternate site energies; the even sites have a lin- 
ear energy e and the odd sites have a nonlinear energy 
x|C n | 2 . Such a system can be described by the Hamil- 
tonian given as 

^=2 E ( C nC* +1 + C*C n+1 ) 
n— — oo 

oo oo 

+ | E l^nl' + f E l^n+ll 4 - (12) 

n— — oo n—~oo 

For x = 0, the system reduces to a periodically modu- 
lated linear system. The translational symmetry of the 
systems are preserved for both x = and x 0. Us- 
ing C n = (j) n exp(-iEt), we obtain the Hamiltonian in 
terms of </>„. It is, however, not possible to find the exact 
ansatz for the localized states in this system, although a 
certain rational ansatz can be considered. For example, 
the onsite peaked solution, intersite peaked or intersite 
dipped solutions are possible. We will consider each of 
these in the subsequent subsections. 

A. Onsite peaked solution 

Since the system has linear and nonlinear sites at reg- 
ular intervals, the solution may have a peak either at 



3 



the linear site or at the nonlinear site. We first consider 
the solution peaked at the linear site, i.e. the solution 
is peaked at the central site of the system. We therefore 
consider the ansatz <j> n — <fo »?'"'■ Using this and the 
normalization condition, we obtain the effective Hamil- 
tonian as 

2n e (I + 77 4 ) X V 4 (l-V 2 ) 2 
cff 1 + /7 2 2 (1 + 77 2 ) 2 2 (1 + r? 2 ) 2 (l - 77 s ) " 

(13) 

The fixed-point equation, dH e ff/dt] = 0, reduces to 

2 (f - 77 8 ) 2 (1 - T? - 677 + 6T7 3 ) - 2x^(1 if f \\ V & ) 
+ 2 X 7 ? 3 (1 - T7 4 ) ((1 - 2 V 2 )(1 T7 8 ) + 277 8 (f - T? 2 )) = 

(14) 

Solving Eq. (14), we found that there are two separate 
regions in the plane of e versus x : onc with three SL 
states another with only one SL state, as shown in Fig. 3. 
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FIG. 4. The phase diagram for the SL states (solutions 
peaked at the nonlinear site) in the plane of e versus \ for the 
system of a one-dimensional chain with linear and nonlinear 
site energies at alternate sites. 
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FIG. 3. The phase diagram for the SL states (solutions 
peaked at the linear site) in the plane of e versus x f° r a 
one-dimensional chain with linear and nonlinear site energies 
at alternate sites. 

The localized solution peaked at the nonlinear site 
may be obtained considering that each site is shifted by 
one lattice site along either positive or negative direc- 
tion, and the resulting fixed-point equation is give as 

2 (l-7y 2 )(l-7 ? 8 ) 2 (l + 7 ? 2 + e7;) 

- XV(1 ~ ?7 4 )(1 ~ V W ) + Xr/(1 - r? 4 )(4vf(l - V 2 ) 

-(l-7 ? 2 ) 2 (l + 7 ? 8 )(l+ 7 7 4 ))= 0. (15) 

The phase diagram for the SL states is presented in 
Fig. 4. We again find that there are two separate re- 
gions, similar to the case peaked at the linear site. Thus, 
at least two onsite peaked SL states (one of them is 
peaked at the linear site and another peaked at the non- 
linear site) may be obtained for any finite values of x 
and e. However, the regions of three SL states peaked 
at the nonlinear site appears from much lower values of 
X and e, compared to the case peaked at the linear site. 



B. intersite peaked and dipped solutions 

For the intersite peaked and dipped solutions we use 
the dimeric ansatz in Eq. (4). After short calculatons, 
we obtain the effective Hamiltonian for the reduced dy- 
namical system as 

H eff = sgn(E)r ] + ^~p 

. e (1 + £V) X (1-V 2 ) 2 (V 4 + P 4 ) (m , 

2 (1 + /? 2 )(1 + V 2 ) 4 (1 + /3 2 ) 2 (1 - ?7 8 ) 1 ' 

with (3 and 77 as defined earlier. Wc first consider the 
symmetric and antisymmetric states, i.e., (3 = ±1. The 
fixed-point equation yields 

4 V ± 

5? = (it# + ^ = / (r/) ' (17) 

where "+" referes to the symmetric states and "-" to 
the antisymmetric states. It should be noted that e has 
no role in the formation of symmetric and antisymmet- 
ric SL states in the system. ,/From Eq. (17) it is clear 
that only one symmetric state appears for any finite 
value of x- However, for antisymmetric states, the sit- 
uation is different. The value of f~{rf) in Eq. (17) ini- 
tially increases as 77 increases and finally decreases as 
77 approaches to 1, leaving a maximum at 77 e [0,1]. 
Thus, there are two critical values of x such that there 
is no anti-symmetric SL state for \ < 18.72, two anti- 
symmetric SL states for 18.72 < x < 32 and one anti- 
symmetric SL state for x > 32, shown as the two vertical 
lines in Fig. 5. 

We next consider the asymmetric states of (3 7^ 1. In 
this case e plays its role, unlike the case of symmetric 
and antisymmetric states. The fixed-point equations are 
given as 

(l + 7 7 2 )(l+/3 2 )(l-7 / 8 ) 2 (l+/3 2 - 2/377) 
-er ? (l-/3 4 )(l-r ? 8 ) 2 - X r?(l-'? 4 ) 2 

X ((T7 4 + /3 4 )(1 + T? 2 + T7 4 - Tf) - 77 2 (1 - T? 8 )) = (18) 
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and 



(l-n 4 )(l+/3 2 )(l-n 8 )(l- 

-^(i + ^xi-^xi-t? 2 ) 

+ X /?(l-n 2 )(l-n 4 )(/3 2 -n 4 ) 



a 2 ) 



0. 



(19) 



Solving Eqs. (18) and (19) we obtain various critical 
lines separating the x — e plane into various regions with 
various number of asymmetric SL states as shown in 
Fig. 5. Complete phase diagram will be obtained by su- 
perposing the Figs. 3, 4, and 5. The figures are shown 
separately because otherwise it will appear very com- 
plicated. The maximum number of states obtained in 
this case is 9. We note that the presence of both linear 
and nonlinear sites increases the number of SL states 
compared to the fully nonlinear case. No SL state is ob- 
tained for x — 0- This is expected because the system 
with x = reduces to a periodically modulated linear 
system which again may be renormalizcd to a perfect 
linear system. 
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FIG. 5. The phase diagram for the SL states due to the 
intersite peaked and dipped solutions in the plane of e versus 
X for the one-dimensional chain with a linear and a nonlinear 
site energies at alternate sites. 



IV. CAYLEY TREE WITH A LINEAR 
IMPURITY AND A NONLINEAR IMPURITY 

Here we consider a Cayley tree with a linear impurity 
and a nonlinear impurity embedded in it. The purpose 
of this work is to observe how the SL states vary accord- 
ing to the geometry of the underlying lattice structure. 
The structure of the Cayley tree for the coordination 
number Z = 3 and the connectivity K — 2 is shown in 
Fig. 6. 




FIG. 6. The structure of Cayley tree with the coordina- 
tion number Z — 3 and the connectivity K — 2. 

We pick a bond and assign the numbers and 1 on 
its two ends. The neighboring sites of the site 1 are 
numbered with an increasing order and those of the site 
are numbered with a decreasing order. Thus, all points 
with the same number are assumed to be in the same 
generation. The number of points in the nth generation 
is K n ~ x if n > 1 and A"H if n < 0. We note that 
all sites in a given generation have the same probability 
amplitude. 

We now consider the motion of a particle on a Cayley 
tree of the connectivity K with two impurities (eo and 
ei) embedded at sites and 1, respectively. In the tight- 
binding formalism with nearest-neighbor hopping only, 
the equations governing the motion of a particle are 



.dC n 
i—— 
dt 

.dC n 
i—— 
dt 

dCt 
i—— 
dt 

.dC 
i—— 

dt 



KC n+1 +C n _ u n>l 
#"C_| n |_i + C_|„| + i, n<0 
KC 2 + C + eiCi 
KCU + C\+e C , 



(20) 



where C n denotes the probability amplitude at any point 
in the nth generation. We assume that the matrix ele- 
ment of the nearest-neighbor hopping is unity and that 
all points in a given generation have the same site en- 
ergy. The normalization condition for the site ampli- 
tudes becomes 



oo 

E ^ h ic„i 2 + -]Tinc„ i 2 = i. 



(21) 



n=l 



Solving the Eq. (20) appears to be nontrivial; how- 
ever, by a simple transformation of variables (i) r = 
VKt, (ii) eo = eo/s/K and ei = c\/\/K, and (iii) 
C n = K-^-V^Cn for n > 1 and C_ H = K^ n ^ 2 C n 
for n < in Eq. (20), we obtain 



.dC n 



G'n+i + C„_i, for n > 1 and n < 0. 
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■ dCi ~ 1 ~ _ ~ 
i-r- = C 2 + -t=Cq + eiCi, 

J-T^ = C_i + -f=C x + e C - 



(22) 



which are precisely the equations for a one-dimensional 
chain, with the normalized probability amplitude C n . 
The normalization condition in Eq. (21) reduces to 
Y^oo I | 2 = 1- Therefore, the motion of a particle on 
a Cayley tree is mapped onto that on a one-dimensional 
chain with a bond defect between the impurity sites. In 
other words, the nearest-neighbor hopping matrix ele- 
ment between the zeroth and the first site is reduced 
from unity to 1/ \[K. It can be shown that the Green's 
function Gofi(E) calculated from Eq. (7) would yield 

Go,o(E = E\[K) for a Cayley tree of the connectivity 
K. Here, we consider a linear impurity at the zeroth site 
and a nonlinear impurity at the first site. Therefore, the 
impurities are defined as ?o = £ and ei = x|Ci| 2 with 
X = X"J~K an d e = e-</K. 

Now, following the same approach as in the case of a 
one-dimensional chain, we obtain all possible symmet- 
ric, antisymmetric and asymmetric SL states for various 
values of the impurity strengths. The full phase diagram 
of SL states in the plane of e versus \ f° r the Cayley tree 
(K = 4) with a linear impurity and a nonlinear impu- 
rity embedded in it is shown in Fig. 7. We note that the 
maximum number of SL states in this case is 7. Thus, a 
linear impurity and a nonlinear impurity play the same 
role as that of the two nonlinear impurities as long as the 
maximum number of SL states is concerned. However, 
the phase diagram is completely different from that for 
a Cayley tree with two nonlinear impurities. The one- 
by-one increment of the number of SL states from one 
region to another is retained. 



V. SUMMARY 

Formation of SL states due to the presence of lin- 
ear and nonlinear impurities is studied. Three different 
cases are considered. First, a one-dimensionl chain with 
a linear impurity and a nonlinear impurity is consid- 
ered, and the maximum number of SL states is found to 
be five. If one considers that the maximum number of 
SL states due to two linear impurities was found to be 
two and that due to two nonlinear impurities was five, 
one may conclude that the presence of one linear im- 
purity and one nonlinear impurity plays the same role 
as that of two nonlinear impurities as long as the maxi- 
mum number of SL states is concered. Variation of the 
energy of SL states as a function of nonlinear strength is 
studied and it is found that when a pair of states of the 
same type appear at the same value of \, on e of them 
goes away from the band edge and another approaches 
the band edge as \ increases. Second, a linear chain 
with linear and nonlinear site energies at alternate sites 
is considered and a rich phase diagrams for SL states 
are obtained. The maximum number of SL states are 
found to be 9, which is larger than that for any sys- 
tem studied so far in this direction. These SL states 
may be called self-localized states because they appear 
even though the translational symmetry in the system 
is preserved. Third, as for the case of a different lattice 
structure, the Cayley tree with a linear impurity and a 
nonlinear impurity is considered and a rich phase dia- 
gram of SL states is obtained. The maximum number 
of SL states is found to be 7. 
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FIG. 7. Phase diagram for SL states in the Cayley tree of 
connectivity K — 2 with a linear and a nonlinear impurity 
embedded in it. 
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